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The groundstate configuration and the eigenmodes of two 
parallel two-dimensional classical atoms are obtained as func- 
tion of the inter-atomic distance (d). The classical particles 
are confined by identical harmonic wells and repel each other 
through a Coulomb potential. As function of d we find several 
structural transitions which are of first or second order. For 
first (second) order transitions the first (second) derivative of 
the energy with respect to d is discontinuous, the radial posi- 
tion of the particles changes discontinuously (continuously) 
and the frequency of the eigenmodes exhibit a jump (one 
mode becomes soft, i.e. its frequency becomes zero). 

PACS numbers:73.20.Dx, 36.40. Ei, 64.90. +b 



In recent years, there has been an increased interest 
in the study of finite two-dimensional (2D) systems of 
charged particles 0. These quantum dots are atomic 
like structures which have interesting optical properties 
and may be of interest for single electron devices. Most of 
the previous works have treated the quantum mechanical 
problem of a small number of electrons. If the number 
of electrons is increased beyond 6-7 the Coulomb correla- 
tion of the electrons has to be treated in an approximate 
way. Therefore, in order to fully understand the physics 
of Coulomb correlations, wc have made a throughout 
study of the classical system in which the particles are 
taken point-like but where no approximation was made 
on the Coulomb correlation. This approach is valid in 
high magnetic fields where the kinetic energy of the elec- 
trons is quenched. In the ground state the charged par- 
ticles are located on rings and a table of Mendeljev was 
constructed in Rcf. (2). These configurations were re- 
covered in a quantum calculation in the limit of large 
magnetic fields ||^. A study of the spectral properties 
of these classical systems such as the energy spectrum, 
the eigenmodes, and the density of states was made in 
Ref. Q|. Such classical confined systems have been ob- 
served in electron dimples on liquid helium in drops 
of colloidal suspensions |^ and in confined dusty parti- 
cles §. 

Here we extend our previous work to the case of clas- 
sical artificial molecules which consist of two classical 
2D atoms which are laterally separated by a distance d. 
Intuitively, we expect interesting behavior as function of 
d which governs the inter atomic interaction. This can 
be seen as follows, for c? = we have just one 2D atom 
(e.g. for 27V = 10, this is a two ring structure), while for 
d — > 00 the system consists of two independent 2D atoms 
with each half of the total number of particles (e.g. two 



atoms each having 5 particles on one ring). This implies 
that as function of d structural transitions (e.g. con- 
figurational changes) have to take place. We found that 
these artificial molecules show a surprisingly complex be- 
havior as function of d. In this letter, we investigate the 
groundstate energy and its derivatives, the groundstate 
configurations and the normal modes, i.e. the eigenfre- 
quencies, of the system. 

We limit ourselves to the system consisting of an even 
number, 2iV, of charged particles which are evenly dis- 
tributed over two layers separated a distance d. In both 
layers, the confinement potential is centered around the 
z-axis and this parabolic confinement, which is taken the 
same for both layers, keeps the system together in the 
xy-plane. We focus our attention on systems described 
by the Hamiltonian (where m is the mass of the parti- 
cles, loq the radial confinement frequency, e the particle 
charge, and e the dielectric constant of the medium the 
particles are moving in) 
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the confinement and interaction energy of the artificial 
atom in layer S = I, II, and 
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(2b) 



the interaction energy between the atoms in the two lay- 
ers. For convenience, we will refer to our charged parti- 
cles as electrons, keeping in mind that they can also be 
ions with charge e and mass m. The Hamiltonian can 
be written in a dimensionless form if we express the co- 
ordinates, energy, and frequency in the following units 
r' = (eVe)i/3a-i/3^^' ^ (eVe)^/^^!/^, = c^o/V2, re- 
spectively, with a = 17110^/2. The numerical values for 
the parameters r' and E' for some typical experimental 
systems were given in Ref. [|| . In dimensionless units the 
Hamiltonian becomes 
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and consequently the groundstate energy is only a func- 
tion of the number of electrons, 27V, and the distance d 
between the layers. 

The numerical method used in the present study to ob- 
tain the groundstate configuration is based on the Monte 
Carlo technique supplemented with the Newton Method 
in order to increase the accuracy of the energy of the 
groundstate configuration. The latter technique is out- 
lined and compared with the Monte Carlo technique in 
Ref. [Q and also yields the eigenfrequencies and the eigen- 
modes of the groundstate configuration. 

Let us consider first the case of a molecule consisting of 
six electrons distributed over two atoms each with three 
electrons. For d — this is a 2D atom with 6 electrons, of 
which we know the groundstate configuration , namely 
(1,5): 5 electrons on a ring and one electron in the center 
of the ring. In the opposite limit, d — s- oo, we have two 
independent 2D atoms, each consisting of three electrons 
for which the groundstate configuration consists of one 
ring containing three electrons |^]. Thus, as function of 
d, we expect a structural transition. 

Fig. shows the energy per electron of the groundstate 
and its first derivative with respect to d. At d = 0.35955 
the first derivative is discontinuous and a structural tran- 
sition takes place. For d < 0.35955 the groundstate 
configuration is (l,2)/(0, 3): the configuration (1,2) in 
one layer (indicated by the open dots in the insets of 
Fig. ^ and (0, 3) in the other layer (indicated by the 
solid dots). As viewed from above we have (1,5) which 
is the configuration of one atom consisting of six elec- 
trons. This implies that the inter-layer correlations are 
sufficiently strong to impose the one atom configuration 
to the 2N — 6 electrons in the molecular structure. For 
d > 0.35955 the configuration is twice (0, 3), which is the 
configuration of two independent atoms each consisting 
of three electrons. 

The frequencies of the normal modes of the ground- 
state configuration are shown in the inset of Fig. |l|. No- 
tice that at the first order transition point the frequen- 
cies exhibit a jump. For d > 0.35955 more modes are 
degenerate in energy and in the limit of d — > oo all eigen- 
frequencies are at least twofold degenerate. The latter is 
a consequence of the fact that the electrons in one layer 
can vibrate in phase and out-of-phase with respect to 
the electrons in the other layer. With decreasing d the 
inter-layer interaction destroys this degeneracy. 

In Fig. [l] also the energies of the metastable states are 
shown. For d — there is just one metastable state with 
configuration (0,6). If d is slightly different from zero, 
there are four metastable states plus the groundstate. 
One obtains these configurations by moving one electron 
from one layer to the other without changing its lateral 
position, and this for each of the two configurations [(1,5) 
and (0,6)]. From Fig. one sees clearly how a configura- 
tion, metastable for d < 0.35955, becomes the stable one 
for d > 0.35955. 

The artificial molecule consisting of two times five 
electrons is far more complex and exhibits several tran- 



sitions some of which are qualitatively different. For 
d — this is a 2D atom with 10 electrons and the 
groundstate configuration is (2,8) Q. For d ^ oo we 
have two independent 2D atoms, both with configura- 
tion (0,5) |2j. Fig. |^(a) shows the eigenfrequencies as a 
function of d and Fig. ||(b) the distance of the electrons 
from the center of the confinement potential. From a first 
glance we have four transition regions, namely around 
d = 0.16, d = 0.68, d = 0.78 and d = 0.90. The spatial 
configurations of the electrons are depicted in the inset 
ofFig.gb). 

Region 1 (around d = 0.16) and region 2 (around 
d — 0.68). The transitions in both regions are quali- 
tatively similar and therefore we limit our discussion to 
d « 0.68. At these transition points there is no abrupt 
change of the configuration, but within a small d-region 
the radii of the electrons change appreciably, but con- 
tinuously (Fig. H(a)). Fig. ^(b) shows that this region is 
delimited by two eigenfrequencies which become zero and 
consequently these transitions are induced by the soften- 
ing of a mode. In Fig. ||(c) the second derivative of the 
energy with respect to d is given. There are two discon- 
tinuities in the second derivative, namely at d = 0.68120 
and d = 0.68215 while the first derivative is continuous. 
In the inset of Fig. ||(c) the spatial configuration is shown 
at d = 0.681 (circles) and at d = 0.683 (triangles). The 
open and closed symbols refer to electrons belonging to 
different layers. Notice that no qualitative changes of 
the configuration occurs at the second order transitions. 
For region 1 we found similar second order transitions at 
d = 0.15775 and d = 0.16175. 

Region 3 (around d = 0.78). The transition occurs at 
d = 0.77605 where the position of the electrons changes 
discontinuously. There is a jump in the eigenfrequencies 
and in the first derivative of the energy. Clearly this is a 
first order transition which is similar to those found for 
the above molecule with two times three electrons. 

Region 4 (around d = 0.90). In this region there 
are two transitions. The first at d = 0.88585, and the 
second at d = 0.90588. Only in the latter the number of 
electrons on the different rings changes. It is a first order 
transition and the eigenfrequencies exhibit a jump. On 
the other hand the transition at d = 0.88585 does not 
correspond to a change in the number of rings, however 
there is a jump in the position of the electrons. It is a 
first order transition with a discontinuity in the frequency 
spectrum in which the smallest non zero eigenfrequency 
decreases for d 0.88585 but stays different from zero. 

For the artificial molecule in which each atom contains 
six electrons we find both first and second order transi- 
tions. With increasing d there are three first order transi- 
tions, namely at d = 0.2965, d = 0.5125 and d = 0.5345, 
followed by one second order transition induced by the 
softening of a mode at d = 0.9259. The latter is different 
from the previous molecule with two times five electrons 
where there are narrow d-regions delimited by two sec- 
ond order transitions. For d < 0.2965 the configuration 
is (2,4)/(l,5), at the first order transitions the position 
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of the particles changes discontinously while at the sec- 
ond order transition the position of the particles changes 
continuously and after it the outer electrons form a per- 
fect circle with the inner electrons sitting in the center, 
resulting in the configuration (1, 5)/(l, 5). 

The groundstate of the artificial atom with seven 
electrons is (1,6), and with eight electrons is (1,7). 
Therefore, we expect that the above new type of tran- 
sition also takes place for the artificial molecule with 
two times seven and two times eight electrons. For the 
molecule in which both atoms contain seven electrons we 
find two narrow d-regions delimited by two second order 
transitions (the first at d — 0.360375 and d = 0.360495, 
the second at d = 0.504625 and d = 0.504865), a first 
order transition a.t d = 0.5425, and indeed the second or- 
der transition at d = 0.9259. Fig. ^ shows the frequency 
spectrum and the position of the particles near the last 
transition. At the second order transition at d = 0.9259 
the system transforms from an imperfect ring structure 
with the middle electron outside the center of the system 
to a perfect ring with one electron in the center of each 
layer (see Fig. ^(b)). 

For the latter molecule another phenomenon occurs. 
Around d « 0.74 the smallest non zero eigenfrequency 
decreases substantially but does not become zero while 
the other eigenfrequencies and the position of the parti- 
cles changes fast but continuously with d. It is clear that 
this is not a first, second or third order transition as is 
made clear in the inset of Fig. ^(b) where we plotted the 
second and third derivative of the energy with respect to 
d. Therefore we refer to this as a continuous transition 
with a lambda-like shape of the derivative d^E/dd^. 

In conclusion, we have presented the results of a nu- 
merical calculation of the groundstate configuration, its 
energy and the spectrum of the normal modes of clas- 
sical double-layer atoms. There is a vertical Coulomb 
coupling between the electrons constituting the molecule. 
These artificial molecules undergo structural transitions 
as function of the distance between the layers. We found 
first and second order transitions. For first order tran- 
sitions the position of the particles changes discontinu- 
ously and the eigenfrequencies exhibit a jump. For sec- 
ond order transitions the position of the particles changes 
continuously and these transitions are induced by the 
softening of an eigenmode. The different transitions are 
summarized in Fig. ^ in a {d, N) diagram for molecules 
containing up to 2A^ = 20 electrons. Although we as- 
sumed that: 1) both atoms of the molecule contain the 
same number of electrons, and 2) both atoms have the 
same confinement potential, a rich variety of structural 
transitions are found. It is expected that if we relax these 
assumptions new configurations are possible. 

The present system constitutes a simple non trivial 
model system exhibiting a rich variety of structural tran- 
sitions. This finite system correctly describes the local 
properties of structural phase transitions. Furthermore, 
we found additionally: 1) transitions which are delimited 
by recurring second order transitions, and 2) a lambda 



type of transition without a divergence in the second 
derivative of the energy. It should be noted that in the 
limit A'^ ^ oo we obtain the system of two parallel clas- 
sical two-dimensional electron gases, each of them form 
a Wigner crystal which was studied in Ref. |^ . 
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FIG. 1. The energy of the groundstate (plotted as dots) 
and its first derivative with respect to the lateral distance d 
between the two atoms constituting the artificial molecule for 
2 A'' = 6. The inset shows the frequency of the normal modes 
of the groundstate. The configurations are schematically rep- 
resented in the figure where solid and open dots correspond 
to electrons in different layers. 

FIG. 2. (a) The eigenfrequencies of the groundstate for the 
artificial molecule with two times five electrons as a function 
of the distance d between the two atoms, (b) The distance of 
the different electrons from the center of the confinement po- 
tential. The three groundstate configurations are also shown 
in the inset of the figure. 
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FIG. 3. (a) The distance of the electrons from the center of 
the confinement potential for the artificial molecule with two 
times five electrons as a function of d. (b) The lowest non zero 
eigenfrcqucncy of the grouridstate. (c) The second derivative 
of the groundstate energy with respect to d. The inset shows 
the configuration at rf = 0.681 (circles) and at d = 0.68.3 (tri- 
angles). Open and solid symbols refer to electrons belonging 
to different atoms. 

FIG. 4. (a) The eigenfrequencies of the groundstate for the 
artificial molecule with two times seven electrons as a function 
of d. (b) The distance of the electrons from the center of 
the confinement potential. In the inset the second and third 
derivative of the energy are shown around d « 0.74. 

FIG. 5. Summary of the different transitions for molecules 
up to 2A'' = 20 electrons. 'Continuous' refers to a transition 
at which the position and the eigenfrequencies change rapidly, 
but continuously, with d and the derivatives of the energy with 
respect to d are continuous. 
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Figure 1 
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